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£3 '. Abstract 

Q 

&JQ' We construct a model of Brownian Motion on a pseudo-Riemannian 

tt . manifold associated with general relativity. There are two aspects of 

. ,— i \ the problem: The first is to define a sequence of stopping times as- 

^" sociated with the Brownian "kicks" or impulses. The second is to 

J^ ■ define the dynamics of the particle along geodesies in between the 

i E~ *. Brownian kicks. When these two aspects are taken together, we can 

associate various distributions with the motion. We will find that 
the statistics of space-time events will obey a temperature dependent 
I/-) , four dimensional Gaussian distribution defined over the quaternions 

CD \ which locally can be identified with Minkowski space. Analogously, 

the statistics of the 4-velocities will obey a kind of Maxwell- Juttner 
distribution. In contrast to previous work, our processes are charac- 
(—■} terized by two independent proper time variables defined with respect 

C*~) ' to the laboratory frame: a discrete one corresponding to the stopping 

times when the impulses take place and a continuous one correspond- 
ing to the geodesic motion in-between impulses. The subsequent dis- 
K^ | tributions are then solutions of partial differential equations which 

contain derivatives with respect to both time variables. 
CS ■ 
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1 Introduction 

Brownian motion is one of the cornerstones of Statistical Mechanics. Ein- 
stein successfully used it to give a rational proof of the existence of atoms 
[lj, and since then our understanding of it has been used as a paradigm 
to model systems in contact with a heat reservoir. It has also been used 
as a stochastic model to represent a variety of different phenomena in such 
diverse fields as physics, chemistry, biology, finance etc. Indeed, its univer- 
sal character rests on it being the simplest model available for describing 
time evolution implied by a combination of random and deterministic fac- 
tors [21 [31 IS]- In the case of Einstein's theory, the deterministic factor is 
given by the Stokes force exerted on pollen grains by a liquid, seen as a 
continuum macroscopic medium; while the random factor represents the 
impulses given to the same pollen grains by the myriad of fluid molecules 
colliding with them. Einstein's ingenuity consisted in understanding that 
water could be seen as acting on pollen grains in two almost antithetical 
ways: as a continuum with its viscosity (systematic component), and as a 
collection of many interacting particles (chance). Clearly, many phenom- 
ena can be interpreted as the result of the cooperation of systematic and 
random events, hence the success of his simple model in combining the two. 
For example, in the case of a tagged molecule of a gas, the free flights be- 
tween collisions with other molecules constitute the systematic part, while 
the chance collisions with other molecules which interrupt the free flights, 
constitute the non-systematic component of the motion causing it to move 
in a random environment [6\\j In principle, this picture applies to all forms 
of dynamics, including special and general relativity. 

One could thus explore the possibility of addressing the relativistic Brow- 
nian motion, as a random process with stationary independent increments, 



Equivalent situations are realized with particles tracing deterministic trajectories in 
regular environments, if correlations decay in time and space, making inapplicable a 
deterministic description. This happens, for instance, in the so-called periodic Lorentz 
gas, consisting of point particles moving in a periodic array of convex (typically circular) 
scatterers, in which position and velocity correlations decay at an exponential rate, [3 [8]. 
Another example is given by polygonal billiards, in which correlations do not decay 
exponentially fast [9] ■ In that case, one observes a different class of phenomena, which 
imply anomalous rather than standard diffusion. 



in which a Brownian particle travels on a geodesic until it is dislodged by 
the interaction with another (moving or standing) object, which shifts it 
onto a new geodesic. 

Seen from this perspective, there are two aspects to consider. One is 
connected with the specific dynamics of the particle along the smooth parts 
of its piecewise smooth trajectory, while the second is connected with the 
random fluctuations that occur as the particle bounces from one smooth 
section of the trajectory to another. In practice, the geodesic motion be- 
tween the interactions obeys the deterministic laws of relativistic dynamics, 
while the collection of impulses assigns random orientations, velocities and 
accelerations according to appropriate statistical laws, although the actual 
accelerations along geodesies will be determined by the dynamics in between 
collisions. 

This perspective will enable us to overcome the historical difficulties as- 
sociated with defining a Brownian motion within the framework of relativity. 
For example, a straightforward extension of the classical theory to the rela- 
tivistic one is known to be frustrated in practice by the difficulty (perhaps 
impossibility) of expressing the dynamics of interacting many body systems 
within the relativistic framework. The choice of reference frame becomes 
problematic especially when accelerations are involved. Furthermore, it is 
hard to formulate the interaction potential because the action and reaction 
principle holds only for contact interactions; and the effects of length con- 
traction and time dilation add to the puzzle [5J [TTJJ [TT]. Also from the point 
of view of stochastic processes and of the diffusion processes associated with 
Brownian motion, with the exception of a non-trivial time-discrete relativis- 
tic Markov model found in [3], certain relativistic generalizations and their 
Gaussian solutions must necessarily be non-Markovian or reduce to singular 
functions (cf. the excellent Review [5J, and references therein). Indeed, this 
has been investigated in great detail by Hakim [12] who defines relativistic 
stochastic processes in \x = A4 A x U 4 where .M 4 is the Minkowski space- 
time and U A is the space of velocity 4-vectors but shows it is not suitable 
for defining a Markov process^ In addition, there is a theorem by Dvoret- 



2 Hakim manages to define relativistic stochastic processes by indexing the stochastic 
events with the set of all regular 7-dimensional space- like hypersurfaces |12| . However, 
that set is only partially ordered, hence not suitable for the Markov property. 



zky, Erdos and Kakutani which states that for Brownian motion the sample 
paths are almost surely nowhere differentiable and one may wonder how 
this fits into the overall framework of relativity and differential manifolds. 

Conceptually, these kind of difficulties can be bypassed, if one is not 
interested in a detailed description of the interaction processes |17[ [H] and 
treats the impulses as a set of random variables occurring at random discrete 
proper timeso 

However, to cast a deterministic process, such as the one described by 
Ref.[T0], in the form of a stochastic process, one cannot index the events by 
their proper time alone. Instead, it is necessary to distinguish the random 
parts from the deterministic parts of the motion by means of stopping times 
and then, in order to be coherent, refer all motions to the same (laboratory) 
frame of reference. Three cases arise. We can have: 

1) A Brownian random walk in which a single time scale is used and 
both the stopping times and the time between jumps are correlated 
by means of a common time. For example, Keller in [21J considers a 
particle moving along the x axis such that during a time interval of 
duration r = 1 from t — i — 1 to t — i, i — 1, 2, . . ., the particle moves 
with velocity v — +1 or v — — 1, each with probability |. 

2) A Brownian random walk in which two independent time scales are 
used, with each stopping time indexed by the time of the Brownian 
kicks as measured in the rest frame of the Brownian particle. In other 
words, the process obeys the strong Markov property and the proper 
time between consecutive stopping times is a random variable, with 
the Brownian flights lying along a piecewise differentiable curve of 
random length each with its own proper time. 



3 For instance, in the low density limit, in which the interaction (potential) energy 
is negligible compared to the kinetic energy, Ref.jTO] describes a relativistic gas as a 
collection of particles which move according to special relativity from collision to collision, 
and treats as classical the "randomly" occurring collisions among particles. In this way 
Ref.[10j provides numerically a dynamical justification of the hypothesis of molecular 
chaos underlying the validity of the relativistic Boltzmann equation and of its equilibrium 
solution known as the Maxwell- Jiittner distribution 1161 1171 . 



3) A continuous time Brownian motion in which two independent proper 
time scales are used with stopping times being chosen such that the 
time between two consecutive stopping times is an arbitrary constant 
(not a random variable) , but with distance and velocity between stop- 
ping times defining a stationary independent Markov process. 

For the purpose of this article, we focus on case (2), leaving case (3) for a 
later date. We propose a new perspective, based on combining the discrete 
stopping times and the continuous time associated with the determinis- 
tic part of the motion. By using the proper time, the same event viewed 
from different frames can be synchronized. Indeed, from the perspective of 
the rest frame (heat bath) of the particle, Brownian impulses recorded at 
Ti, T2, . . . Tj . . . can be transmitted to the laboratory frame. For example, 
the numbers .5, .23, .678, .45, ... could represent the times between succe- 
sive impulses as recorded in the rest frame, which for a large sample can be 
used in principle to calculate the variance of n — Tj_i. Since each r, — Tj_i 
is a proper time (recall the local time of the rest frame is equivalent to 
the proper time) it remains invariant in general relativity under coordinate 
transformations, which also includes the laboratory frame. This is the only 
information transmitted from frame to frame. Moreover, if we assume that 
the heat bath is held at a constant temperature, then it is also reasonable to 
assume using the strong Markov property that the set {tj — Tj_i}^ 1 defines 
a set of stationary independent increments. This statement is independent 
of reference frames. 

It is also true that from the laboratory frame's perspective we do not 
know the local coordinates (t, x, y, z) of the Brownian particle and that 
many options are possible. However, we can extract enough information 
from the invariance of the {tj — Tj_i} and the Central Limit Theorem to 
specify the probability of finding the particle in any given region defined 
with respect to a reference frame. Our density function will be Lorentz 
invariant (see below) and consequently allows us to calculate probabilities 
in any reference frame. Indeed, that is why we are formulating the theory 
from the statistical perspective and not deterministic mechanics. 

To conclude, the novelty of this article will be to show that a relativis- 
tic Brownian motion can be described when motion between the discrete 
random stopping times lies along a geodesic path. In doing so, we obtain 



a covariant description of the relativistic Brownina motion, analogously to 
Ref.jTS], and we also circumvent the difficulty in dealing with nowhere dif- 
ferentiable paths associated with case (3). 

For the model developed here, we show that in the case of general or 
special relativity, Brownian motion can exist along piecewise geodesic curves 
on a pseudo-Riemannian manifold. 

2 Working with geodesies 

The key to the above mentioned developments will be found in working with 
motion along geodesies. Beginning with the metric tensor 

ds 2 = Qijdx 1 dx^ (1) 

this can be re-written as 

dx 1 
ds = QiA——dx 3 . 
ds 

However, from the definition of a geodesic, we can erect a tetrad {e^} such 

that 

dt dx dy dz 

ds = c—dr —dx — ay —dz (2) 

dr ds ds ds 

along a geodesic. Sometimes, we will prefer to write this in the form 

eds = Hdt — pidx 1 — p2dx 2 — p^dx 3 (3) 

where p is denoted by —H, c is the speed of light in a vacuum and pi = 
c(dx l /ds). Note that since ds is an exact differential s has the form of a 
Hamilton- Jacobi function. In particular, if we change the parametrization 
from t to t* then equation ([3]) can be rewritten as 

eds* = H*dt - pidx 1 - p* 2 dx 2 - p* 3 dx 3 . (4) 

Equations ([3]) and @ can be interpreted in two ways. From one perspec- 
tive they represent two different geodesies parameterized with respect to a 
common time parameter; from another they can be seen as representing the 
same curve with respect to different tetrads. The important thing is that 
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in either interpretation we can pass from one representation to the other 
by means of the expression ds = -^ds*. For our purposes, it is sufficient to 
erect a tetrad in the rest frame of the particle and define a time parameter 
r with respect to a standard clock in this frame. 

In the case of a piecewise geodesic ds = J^ dsi where each dsi defines a 
line segment on a geodesic with respect to the laboratory tetrad, each Sj can 
be parameterized in terms of S\, say, which in turn can be parameterized by 
r. In other words, each component Si is a Hamilton- Jacobi function which 
can be synchronized with the laboratory frame by means of the expression 

dsi dsi ds\ 
dr ds\ dr 

Note -jj: does not have to be a constant. In particular, if r is chosen as a 
standard time parameter then all geodesies can be parameterized with re- 
spect to the same parameter by means of this relationship. This is equivalent 
to defining a universal time parameter [18], first introduced by Stuekelberg 
in 1941, [19J, and further developed by Horwitz et al. [20J. 

With a universal laboratory time established, we now derive two other 
equations from the metric which will be useful later on for describing spe- 
cific realizations of Brownian processes. Equation ([I]) expressed in tetrad 
notation becomes 

ds 2 = c 2 dt 2 — dx 2 — dy 2 — dz 2 . 

Letting s = dsjdr and s = d 2 s/dr 2 this can be rewritten with respect to 
the time parameter r as 

sds = c 2 idt — xdx — ydy — zdz . (5) 

In the case where the standard time r is a non-affinej parameter of the local 
time t, the second derivative does not vanish and we can write 



') 



or equivalently 



ss = c 2 ii — xx — yy — z'z (6) 

sds = c 2 idt — xdx — ydy — 'zdz (7) 



For an afRne parameter there is no acceleration which means s = 0. 



or 

sds = c 2 idt — xdx — ydy — zdz, (8) 

where equation (8) follows from (7) by observing that for any function / 

A A 2 f d df ■ ■ 

dr 1 dr dr 

Note that since sir) is a scalar function, sds = ds 2 is both Lorentz invariant 
by construction and also an exact differential with respect to r. 

3 Brownian Motion 

To investigate a possible notion of relativistic Brownian motion, let us return 
to equation (jSJ) 

sds = c 2 idt — xdx — ydy — zdz , 

defined with respect to a universal time parameter r at any point along a 
curve. For what follows, we will restrict the motion to piecewise geodesic 
curves with Varfc — Ti_i) > 0; in other words, to those curves for which the 
tangent is defined at every point along a (smooth) piece of the trajectory 
geodesic in M. in between two instantaneous random events occurring at 
discrete times r»_i and n, which deviate it from its free motion_| such that 



re-l 

s(r) = 



^ / dsi+ ds n . 



To put this into the context of probability theory, we now need to define 
random variables on the sample space Q = {s} of timelike geodesies on M.. 
Let S, S and S be three random variables on Q defined by 



ds » d 2 s 

Tr ^ S{S) = d? 



S(s) = s, S(s) = — and S(s) = -y-^ with r G [0, oo) 



On occasion, we will also write S = (T, X, Y, Z) or S = X a , a e {0, 1, 2, 3} 
where each X a is a projected cosine of S with respect to a tetrad. In 



5 because of collisions with other particles or with various kinds of obstacles [HI El [9] . 



addition, we can define a random differential dS by dS(s(r)) = ds and a 
random variable 

SdS = c 2 fdT - XdX - YdY - ZdZ = X a dX a , with r G [0, oo) , (9) 

which for each event s can be associated with the metric (151) . Also, in the 
event that motion along the geodesic is not at a constant speed with respect 
to the standard clock, a random variable associated with acceleration can 
be identified with equation (7) for each realization of the motion. It is given 
by 

SdS = c 2 fdT - XdX - YdY - ZdZ = X a dX a (10) 

We now introduce a set of stopping time {r^} and a set of random vari- 
ables Si, Si, Si, indexed by the stopping times, and defined over the sample 
space of timelike geodesies {s(r)} such that for all r G [tj_i, Ti) 

ds •• d s 

Si(s) = sin- Ti_i), Si(s) = -i-(Ti) and Si(s) = ~p^( T i)- 

In particular, in the case of a piecewise smooth curve a(s) = |Jo"j(s) 
where <Jj(s) is a line segment lying on a geodesic with differential dsi, indexed 
by the stopping times, the overall length of a trajectory of a particle under- 
going Brownian motion during the time interval [0, r) where r G [r n _i, r n ), 
can be re-expressed as: 



"- 1 rS Ti r-S T n-1 .s T . ,S T 

s{r) = z~2 I ds i+ ds n = ^2 dS i (a i (s))+ / dS n (cr n (s)). 
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It follows that for each realization of the processes {Ti}°l ± defined with 
respect to the standard clock in the rest frame of the laboratory, {sj}^ 1; 
{sj}^ 1 and {s,}^ 1 define Markov random walks with respect to the (well 
ordered) index set of stopping times rj. In all we have at least five different 
random variables indexed by the stopping times: the difference between 
stopping times (xi = T i ~~ T i-i)i the length traveled between stopping times 
(Si), the initial speed \V®\) at each stopping time, the initial acceleration 



magnitude (|-F"|) at each stopping time and the unit velocity (direction of 
the trajectory) ( ppre a ) with respect to the laboratory inertial frame at 

each stopping timejj 

Suppose now that the heat bath is characterized by a unique constant pa- 
rameter, the temperature B, as in the classical case. Then all five families 
of random variables can be considered to be stationary independent pro- 
cesses and the specific form of the resulting Markov processes will depend 
on the distributions of these random variables. Also, in this case {Si}^, 
{Si}^ and {Si}^ will have infinitely divisible distributions which can be 
associated with Levy processes. 

The details of the relevant statistics depend on the properties of the 
environment in which the process takes place, commonly known as the heat 
bath. As usual in Brownian motion theory, we assume that the only relevant 
parameter of the heat bath is its temperature 6. Consequently, for constant 
temperature 9,5* = YH=i Si, S = XT=i ^i and 5 = Y^l=i &i converge in 
distribution to stable processes. Indeed, this 0-dependence will influence 
the time increments \i an d the measure of a universal unit of time. In 
what follows, we assume (as suggested by molecular dynamics studies such 
as [TU1 [TT])that the underlying thermal bath state is characterized by the 
isotropy of test particle trajectories associated with the independent and 
identically distributed (with 0-dependent distribution) time-like increments 
of stopping times \i- 

With notation clarified, we now construct a statistical model for rela- 
tivistic Brownian motion by examining two cases[] 



6 In the geodesic parts of the motion the direction of velocity equals the direction of 
acceleration. 

7 As it turns out, in any inertial frame an infinite number of Markov processes (one 
for each k) can be defined along a geodesic by 



d k S t 
dr k 



d Si d Si— i . l-ii *st 

are stationary and independent, i£N 



dr k dr k 



provided the derivatives are not everywhere equal to zero. 
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3.1 Relativistic Gaussian Distribution: 

Consider the family of stopping times {tj : i 6 N}, assuming that the 
set of random variables {\i = t% — i~j-i} are independent and identically 
distributed. The fundamental metric is 

ds 2 = c 2 dt 2 — dx 2 — dy 2 — dz 2 

which can be differentiated with respect to the parameter r yielding Eq.([5]). 
This metric can be associated with the quaternions by means of the inner 
product 

(<?i, 92> = ^(<?i<?2 + q{q* 2 + <?2<?i + q* 2 ql) = x y - £i2/i - x 2 y 2 - x 3 y 3 

which is equivalent to the Minkowski inner product, where q\ and q 2 are 
quaternions and q{,q 2 their conjugates. Moreover, since each quaternion is 
also a four vector on Minkowski space, we can define a Lorentz transforma- 
tion q' = Aq such that 

(q[,q 2 ) = x' y' - x[y[ - x' 2 y' 2 - x' 3 y' 3 

= x y — X\y\ — x 2 y 2 — x^y^ by Lorentz transformation 

= (<?1,<?2)- 

Indeed, written in this way the quaternions together with this inner 
product can be identified with Minkowski space. This now enables us to 
associate a random variable Si = S(xi) and random vectors 5* = Tl + Xj + 
Yk + Zl with each stopping time r, such that 

{(T, X, Y, Z) e Q\dS 2 = c 2 dT 2 - dX 2 - dY 2 - dZ 2 } 

where 1, j, k, I are the standard basis for Q, the set of quaternions. If between 
stopping times a particle is moving with velocity Vj(t) = (ct,x,y,z) G Q, 
along a geodesic then the distance traveled componentwise between stopping 
times (tj_i,Tj) is given by 

x( Xl ) a = I l v a dr > 0, (12) 

J Ti-\ 
11 



which, because of the randomly assigned initial velocities at each Tj, can be 
rewritten as a stochastic random vector (quaternion) 

Xf = f VU{r)dr, (13) 

where Xf = X a (xi). It follows by isotropy that for each a G {1,2,3}, 
M{Xf) = lE(X?e a ) = and it also follows from isotropy and (JT2D that 
Var(X}) = Var(Xf) = Var(Xf) > exists. In particular, if V£i(s) is 
constant along s then 

Xt = VUXi, (14) 

and 

Var(Xf) = Varivi-jVarfa - 7i_i) = a»a 2 ( X ) (15) 

exists. 

However, it is not immediately clear how we should interpret the time 
variable. First note that Xi = T i ~ r «-i = — fa-i — Tj). Seen from the 
perspective of the event Tj_i, x« equals the time between the z — 1 kick and 
the future kick at r». Therefore, lE(\xi\) = fit is the mean future time to 
the next event. However, seen from the perspective of r», Xi equals the time 
between the % kick and the previous kick at Ti_i. Therefore 2E7(— |x,|) = — /x t 
represents the mean past time to the previous kick. Consequently it is 
meaningful to ask what is the probability that the next kick will take place 
within t seconds and it is also meaningful to ask what is the probability 
that the previous kick took place within the past t seconds. Taking the two 
statements together, we can maintain the arrow of time going from past to 
future, while meaningfully extending the range of x to (— oo, oo). For this 
extended range of Xi it again follows by symmetry that lE(xi) = 0. 

In particular for the random variable vector 

n n 

s* n = J2 x ? = Y, s t ( 16 ) 

composed of a sum of i.i.d. variables X" — X" _ we can associate a proba- 
bility measure over the quaternions given by 

dP(n,x a ) = P(S^ e (ct,ct+cdt)x(jx,jx+jdx))x(ky,ky+kdy)x(lz,lz+ldz)). 

12 



By the central limit theorem S* is asymptotic for large n to the (quaternion) 
Gaussian distribution in Minkowski space, with variance growing linearly 
with n, given by 

1 f 1 / c 2 t 2 x 2 y 2 z 2 

f s (tl+jx+ky+lz;n) ~ exp 



47T 2 n 2 0" x crxcrycr2 [ 2n \ <7 2 0"^ <7y <7§ 

(17) 
which holds for c 2 t 2 — x 2 — y 2 — z 2 = 0(n) or less. 
As noted above, for the quaternion 

1 / ct . x y z 

q = — 1 j k l- 



2n V a x a x &y ®z 

•s 2 = (?) q) is invariant. This can also be expressed in tensor notation by 
x l = a % Ax\ where A = a* is a Lorentz transformation such that: 

S — — z^ij^ijiJii X — — CLj CL ■ 2-j<ijCLj/ZC CL-/CC — — 2~i*ijOC OC 

It can also be written more succinctly as 

s 2 = {q') j Z'q' = (AqYAEA^Aq = g f Sg, 
where £ is the covariant matrix of the vector in Minkowski space. 

3.2 Diffusion Equation 

Before writing down the diffusion equation for the above process, it is im- 
portant to note that there are two time parameters expressed in Eq. ffT7|) 
and both are independent of each other. The continuous \ parameter is as- 
sociated with the systematic geodesic motion between two stopping times, 
the other with the randomly occurring discrete set of stopping times {t^}. 
The first is directly related to the proper times between stopping times and 
is used to index the deterministic part of the process, while the index n 
associated with the stopping times is used to index the relativistic space- 
time stochastic processes. Once these distinctions are kept in mind the 
usual paradox associated with interpreting random walks, whereby there 
is a finite probability that the Brownian particle may be in superluminous 
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regions, no longer exists. Indeed, Eq. (j!7p was constructed with the assump- 
tion that all deterministic motion is along timelike geodesies. Consequently 
the motion can never be in a superluminous region. 

With this distinction made, one observes that a relativistic Brownian 
particle obeys the following kind of "diffusion" equation: 

al^F _ a%&F _ a^d^F_ _ aJd^F_ _ d_F_ = 

2c 2 dt 2 2 dx 2 2 dy 2 2 dz 2 dn l ' 

This formally resembles the telegraph equation on Minkowski space as found 
for example in Ref.[2T]. However, it is radically different in that the first 
and second order time derivatives of (fT8|) are taken with respect to different 
time variables. Moreover, there is no bound on the occurrences of stopping 
times. The time and the distance traveled between stopping times may be 
arbitrarily large, as long as they are consistent with speeds smaller than 
c. This generalizes the approach of e.g. Ref.[21j, in which the stopping 
times are deterministically generated and equal n at the n-th collision, for 
a one dimensional random walk that increases or decreases incrementally 
by one unit during each time interval. In other words, the case discussed 
by Keller is a special instance of our time-like condition, in which our two 
time variables become a unique variable, leading to the standard telegraph 
equation for the propagation of the particle. 

3.3 Maxwell- J iittner Distribution: 

If instead of Eq.([S]) we begin with the metric Eq.([H]): 

sds = c 2 idt — xdx — ydy — zdz 
we can now repeat the above theory by replacing Eq.( fT3l) with the equation 

X? = V{r t ) - V(rr) + r F^dr (19) 

where V{ji) — V(rf) defines the random discontinuity in velocities that oc- 
cur at each stopping time and F°_ ± represents the four-acceleration along the 
geodesic segment from time 7"i_i to time Tj. For example F^_ 1 could be the 
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acceleration of a Brownian particle performing a simple harmonic motion 
betweens kicks. The above can now be repeated to generate a distribu- 
tion which formally looks like the Maxwell- Jiittner distribution of special 
relativisitc gases: 



1 / 1 ,c t x y z 



f v (lt+jx + ky + lz) ~ exp - — (— r-^V- — ) ). (20) 

(27m) 2 |S| V In a z f a^ a^ a\ J 

but because of the two time parameters t and n, it is radically different from 
the Maxwell-Jiittner distribution, and it applies to accelerating particles. In 
any event, /„ is clearly covariant since |E| = OfO^GyO^ is invariant. 

An interesting special case occurs when there are no accelerations along 
the geodesic and velocity is a constant. In this case, equation (19) reduces 
to 

AX? = V(n) - V(r-) (21) 

and the random variable S* = Y^h=\ V? defines a sum of four independent 
and stationary random vectors over the quaternions and converges to the 
distribution given in equation (20). Clearly the same argument can be 
extended to any random variable such as acceleration that undergoes a 
change from r~ and Tj. 

The above derivation suggests that at the core of the molecular chaos 
hypothesis in relativity are the stationary independent random velocities 
(and accelerations) produced at stationary and independent random stop- 
ping times. This would need to be investigated more in depth, since it 
departs sensibly from the standard wisdom on the Brownian motion. 



4 Concluding remarks 

1. It is important to emphasize that in our calculation in order to have a 
probability distribution, integration is carried out over the quaternions 
and not R 4 . For example 



t\ rib 



»2 Jl 



I />ti no _ t z , x * 

P(T e(h,t 2 ),X e(a,b)) = / / e ^ ^dtd{ix) 

27 rnJ tn J ia 



1 



to Jia 
b 
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to Jo, 



e ^^r 2na x dtdx 



This is equivalent to have constructed a Markov chain on a Minkowski 
manifold subjected to the conditions of special relativity. 

2. In the above formulation we have extended this to pseudo-Riemannian 
manifolds by defining the geodesic with respect to a fixed tetrad, which 
means that locally in can be identified with Minkowski space. 

3. In both the Brownian motion and the Maxwell Jiittner distribution, 
the above theory can be extended to non-geodesic smooth curves W(s) 
parameterized by the geodesic curve length s by means of the relation 
dW = ^ds. 

as 

4. Clearly f s (li + jx + ky + Iz) is invariant under the (homogeneous) 
Lorentz group, while f v (H + jx + ky + Iz) is invariant under Poincare 
group. 

5. The distribution f s and f v can be factorized into separate time and 
space component products. For example 

/(I + jx + ky + Iz) = f n (t)fx(jx + ky + Iz), 

where f^_(jx + ky + Iz) is Gaussian. 

The above, leads to a characterization of the equilibrium state in the 
general relativity framework which, along the lines of classical statistical me- 
chanics, requires almost no information about the details of the interactions 
among the objects of interest. Nevertheless, these details are important and 
implicit in the above formulation. In the first place they are the underlying 
cause of the unpredictable motions associated with the random interaction 
times and are constitutive of the phenomenon. In the second place, the 
statistics of these interactions depend on the temperature G. Knowledge 
of these details is implicitly contained in the relation between the stopping 
times and the laboratory frame time, and should be made explicit to fully 
characterize the relativistic Brownian motion. 
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